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Abstract 



We studied a time history of X-ray spectral states of a black-hole candidate, 4U 
1630-47, utilizing data from a number of monitoring observations with the Rossi X- 



Ray Timing Explorer over 1996-2004. These observations covered five outbursts of 
4U 1630-47, and most of the data recorded typical features of the high/soft states. 



We found that the spectra in the high/soft states can be further classified into three 
states. The first spectral state is explained by a concept of the standard accretion disk 
Oh" picture. The second state appears in the so-called very high state, where a dominant 

hard component is seen and the disk radius apparently becomes too small. These 
"55 ! phenomena are explained by the effect of inverse Compton scattering of disk photons, 

as shown by Kubota, Makishima, & Ebisawa (2001, ApJ, 560, L147) for GRO J1655- 
40. The third state is characterized in such a way that the disk luminosity varies in 
proportion to T ; 2 , rather than T^, where T in is the inner-disk temperature. This state 
is suggested to be an optically-thick and advection-dominated slim disk, as suggested 
by Kubota & Makishima (2004, ApJ, 601, 428) for XTE J1550-564. The second and 
third states appear, with good reproducibility, when T in and the total X-ray luminosity 
are higher than 1.2 keV and ~ 2.5 x 10 38 (£>/10 kpc) 2 cos6>/ (l/\/3j erg s _1 , 
respectively, where D is the distance to the object and 9 is the inclination angle 
to the disk. The present results suggest that these three spectral states commonly 
appear among black-hole binaries under high accretion rates. 

Key words: accretion, accretion disks — black hole physics — stars: individual 
(4U 1630-47) — X-rays: stars 
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1. Introduction 



The physics of accreting black holes is naturally the same among stellar-mass black 
holes, ultraluminous X-ray sources (ULXs), and active galactic nuclei (AGN), and information 
concerning accretion disks enables us to understand the spacetime geometry around the central 
black holes. Among many classes of black holes, galactic black hole binaries (hereafter BHBs) 
are the best for detailed spectral and timing studies, and often the observational results of BHBs 
are considered as a basis to understand other classes of black holes. Therefore, it is essentially 
important to precisely understand the accretion-disk structure of BHBs from observed spectra. 

In classical understanding, the X-ray spectrum of BHBs exhibits two distinct features 
of high/soft and low/hard states (e.g. Tanaka, Lewin 1995). In the high/soft state, the 
spectrum consists of a dominant soft component accompanied by a powerlaw component with 
a photon index of T = 2-2.5 (e.g., Tanaka, Lewin 1995). The dominant soft component 
is believed to be a thermal emission from a geometrically thin and optically thick standard 
accretion disk (Shakura, Sunyaev 1973) extending down to the last Keplerian orbit around the 
central black hole. This is the standard view of the high/soft state BHBs. In fact, this soft 
component is well reproduced by a multicolor disk model (MCD model or diskbb in the xspec 
package; Mitsuda et al. 1984) that approximates the integrated emission from the standard 
disk. The MCD model has two parameters: an inner disk temperature T m and a normalization 
Kmcd — r in/(-D/10 kpc) 2 cos# with an apparent disk inner radius r in , a disk inclination angle 
6, and a source distance D. Correcting r- m with the inner boundary condition (Kubota et 
al. 1998) and color-to-effective temperature, we can estimate the true inner radius, R[ n . Past 
observations of many BHBs revealed that Ri n is kept remarkably constant even though the 
X-ray luminosity varies considerably (Makishima et al. 1986; Takizawa 1991; Ebisawa et al. 
1993, 1994). Moreover, through studies of many BHBs, including Cyg X-l (Dotani et al. 1997), 
the inner radius R m is found to be consistent with the last stable orbit, R so = 3R g = QGM/c 2 , 
where R g is the Schwarzschild radius, G is the gravitational constant, c is the light velocity, 
and M is the mass of the non-spinning black hole. Therefore, it is strongly supported that the 
standard view is realized in the high/soft state BHBs. 

However, this simple, but physically well-understood, picture was found to be sometimes 
broken when the disk luminosity Ldisk reaches a certain critical luminosity. The biggest problem 
is that the estimated values of i?; n do not remain constant when the powerlaw contribution be- 
comes significant for the total luminosity. Based on the Rossi X-Ray Timing Explorer (RXTE) 
observations of the galactic BHBs GRO J1655 — 40 and XTE J1550 — 564, Kubota, Makishima, 
Ebisawa (2001), Kubota, Makishima (2004b) found that this anomaly is observed when the 
disk luminosity exceeds a certain critical value (or a critical temperature around 1.0-1.2 keV), 
while the standard view is realized below such a luminosity. They called this branch an anoma- 
lous regime to distinguish from the standard high/soft state (standard regime). Together with 
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the anomalies in i?j n , the anomalous regime can be characterized by a large intensity variation 
associated with quasi-periodic oscillations (QPOs) and the strong hard emission described by 
steeper powerlaw of T = 2.5-3. The characteristics of the anomalous regime are naturally the 
same as those of the very high state in the literatures (e.g., Miyamoto et al. 1993; McClintock, 
Remillard 2003; van der Klis 2004). Kubota, Makishima, Ebisawa (2001), Kubota, Makishima 
(2004b) found that the strong hard emission is caused by inverse Compton scattering due to 
high-energy electrons that may reside around the disk, and for the first time found quantita- 
tively that the i?i n returns to the same value as seen in the standard regime, if the Compton 
effect is taken into account. This scenario is consistent with those of QPOs, which are usually 
thought to be related to the corona surrounding the disk. 

Kubota, Makishima, Ebisawa (2001), Kubota, Makishima (2004b) also found that 
another unusual state appeared in the brightest period of XTE J1550-564 and the beginning of 
the outburst of GRO J1655 — 40. The X-ray spectrum is apparently expressed by the dominant 
soft component and very weak powerlaw tail like in the case of the standard high/soft state. 
However, Ldisk is not proportional to T^, but to (or r in oc T;~ ) under the MCD plus powerlaw 
model, and the spectral shape of the soft component is somewhat distorted from that of the 
standard disk. This regime is called an apparently standard regime by Kubota, Makishima 
(2004b). They suggested that the apparently standard regime shows a sign of a slim disk, in 
which optically thick advective cooling becomes significant (Abramowicz et al. 1988). 

These observed three spectral regimes seem to be consistent with the theoretically pre- 
dicted S'-shape sequence for the optically thick accretion-disk solutions. This correspondence is 
simple, and it is sometimes referred to understand the spectral behavior of ULXs (e.g., Kubota, 
Done, Makishima 2002; Mizuno et al. 2001). However, such three regimes have ever been con- 
firmed for only two Galactic black holes, GRO J1655-40 and XTE J1550-564, and it is not 
yet clear whether these regimes are general. Moreover, we have not known whether the critical 
luminosity is always the same in recurrent BHBs that repeat outbursts. In this paper, we hence 
present spectral analyses of a recurrent black hole binary, 4U 1630-47. 

This object is one of the famous historical BHBs (e.g., Tanaka, Lewin 1995), and is 
known to repeat X-ray outbursts with a period of about 600-650 days (Jones et al. 1976; 
Parmar et al. 1995). No optical counterpart is known for 4U 1630-47, probably due to its 
high reddening and crowded field (Parmar et al. 1986). A number of authors have shown 
the observational results of this source based on X-ray observations with Ginga (Parmar et 
al. 1997), EXOSAT (Parmar et al. 1986; Kuulkers et al. 1997), BATSE/CGRO (Bloser et 
al. 1996), Beppo-SAX (Oosterbroek et al. 1998), and so on. We here analyzed 322 RXTE 
data of this source during 1996-2004 observations, consisting of 5 outbursts. RXTE data of 
the 1996 and 1998 outbursts have already been analyzed by other authors, including Tomsick, 
Lapshov, Kaaret (1998), Kuulkers et al. (1998), Cui, Chen, Zhang (2000), Hjellming et al. 
(1999), Dieters et al. (2000), Trudolyubov, Borozdin, Priedhorsky (2001), and Kalemci et al. 
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(2004). The anomaly of R m (e.g., Oosterbroek et al. 1998), together with existence of QPOs 
(e.g., Tomsick, Kaaret 2000) during the 1998 outburst has been reported. 

The aim of the present work was (1) to test whether the anomalous regime and apparently 
standard regime are found in the complete RXTE data set of this source, (2) to determine the 
critical luminosity and to test whether the critical luminosity is the same among all outbursts 
of this single source, (3) to confirm the scenario of the inverse Compton scattering in the 
anomalous regime (i.e., to test whether the scenario can solve the problem of variable -Rm), and 
(4) to test the slim disk scenario if the apparently standard regime appears, and then we can 
discuss the generality of the scenario suggested by Kubota, Makishima, Ebisawa (2001) and 
Kubota, Makishima (2004b). This is very important to understand the spectral evolution of 
accreting black holes. The black hole mass, distance D, and inclination 9 of 4U 1630-47 are 
not determined; in this paper we assume D = 10 kpc and cos# = l/y/3, following the previous 
works, and do not treat the data in the low state. 

Before presenting the analyses, we briefly compare the spectral regimes applied in this 
paper to a more seminal definition of the spectral and timing behavior used in other literature in 
section 2. In section 3, we show the observation and data-reduction criteria, and in subsection 
4.1, we characterize the observed spectral states based on the canonical MCD plus powerlaw 
model. In subsection 4.2 and subsection 4.3, we consider the effects of the inverse Compton 
effect and the slim disk, respectively. In section 5, we give a discussion on the results from the 
view point of the disk structure. 

2. Brief Notes On Seminal Definition of Spectral/Timing States 

The spectral/timing states of BHBs are sometimes complicated. Therefore, we here 
briefly show the correspondence of the spectral states that we used in the present work to the 
states defined in other literature. A summary is given in table 1. Beyond the two classical low 
and high states, several subdivisions are defined in the literature: the very high state and/or 
the intermediate state, the disk dominant high/soft state, the low/hard, and the off/quiescent 
state (e.g., McClintock, Remillard 2003). 

The off/quiescent state is in an extraordinarily faint period (Lx — 10 30 - 5-33 - 5 erg s _1 ) and 
generally understood by a picture of advection-dominated accretion flow (ADAF: Narayan, Yi 
1994, 1995; Narayan et al. 1996). The low/hard state is characterized by a single powerlaw 
spectrum of ~ 1.7 with an exponential cutoff around several tens of keV; the accretion is 
thought to still be dominated by advection. Large intensity variation of a strong flat-top noise 
accompanied by type-C (or C) QPOs with a variable centroid frequency (Remillard et al. 
2002) is usually observed in this state (McClintock, Remillard 2003). 

The disk-dominant high/soft state is just the same as the standard regime, where the 
standard disk is realized. The intensity variation is very weak and QPOs are seldom observed, 
except on rare occasion of type-A QPOs associated with a weak red noise (e.g., Remillard et 
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Table 1. Correspondence of the definition of the states in this paper to those in other papers. 



this paper 



classification 



MR03* QPO type 



t 



Remarks * 



(off/quiescent) off/quiescent 
(low/ hard) low/ hard 

(very high state very high (intermediate) 
with strong 



hard emission) 
anomalous 
standard 
apparently stan- 
dard 



very high (intermediate) 
high/soft 
high/soft 



hard 
hard 
SPL 



SPL 
TD 
TD 



unknown 
C or C 
C or C 



B 

A or none 
none 



thermal IC 
thermal IC 
thermal IC 



thermal IC 
disk dominant + PL 
disk dominant + weak PL 



* State definition by McClintock and Remillard (2003). SPL and TD represents a steep powerlaw state 
and a thermal dominant state, respectively. 

t QPO type is referred to Remillard ct al. (2002). Note that majority of QPOs are found as type B. 
■f Characteristics of X-ray spectrum. IC and PL represent an inverse Compton scattering model and a 
powerlaw model, respectively 

al. 2002; Nespoli et al. 2003; Casella et al. 2004; Homan et al. 2001). Based on the 
spectral shape and the luminosity, the apparently standard regime is also classified into the 
disk-dominant high/soft state, even though we are suggesting that the standard disk is not 
justified in this regime. 

Definitions of the very high state and/or the intermediate state are somehow controver- 
sial. In the original paper by Miyamoto et al. (1993), the very high state was first defined 
as a high luminosity state with a large intensity variation, in contradiction to the standard 
high/soft state. After their pioneering work, spectral studies revealed that the very high state 
is associated with a very strong hard component with a steeper photon index of V — 2.4-3.0 
than the hard tail in the standard high/soft state (McClintock, Remillard 2003 and references 
therein). Based on this strong and steep hard component, McClintock, Remillard (2003) called 
the very high state the steep powerlaw state. Detailed timing studies found that, in the very 
high state, QPOs were generally observed like in the case of the low/hard state (van der Klis 
1994, 2004 and references therein). In the case of XTE J1550-564, type-B QPOs are observed 
when the hard emission is less dominant (but still much stronger than the standard high/soft 
state), while type-C (or C) QPOs are observed when the hard component is the most signif- 
icant in the total emission (e.g., Remillard et al. 2002; Kubota, Makishima 2004b). Some 
authors suggest that the very high state is itself the intermediate state of transition between 
the low/hard state and the disk-dominant high/soft state (e.g., Rutledge et al. 1999; Homan 
et al. 2001; van der Klis 1994), while other authors suggest that the very high state is a bona 
fide state of BHBs, and note that nothing in their definitions constrains the order in which 
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state transitions should occur (e.g., McClintock, Remillard 2003). Recently, Belloni et al. 
(2005) classified the very high state and/or intermediate state into two regimes along the flare 
evolution. As previously described, the anomalous regime is basically the same as the very high 
state, but strictly speaking, corresponds to the very high state with type-B QPOs and weaker 
hard emission (Kubota, Done 2004a) or the soft intermediate state (Belloni et al. 2005). The 
very high state with much stronger hard emission accompanied by type-C QPOs was found just 
after the low/hard state, would correspond to the hard intermediate state (Belloni et al. 2005), 
and the optically thick disk is somehow truncated at > R so (Kubota, Done 2004a). Thus, this 
type of the very high state may be associated with the state transition. 

Thus, the states of bright BHBs are still being discussed from the view point of both 
spectral and timing phenomena. On the other hand, three regimes treated in the present work 
are based on the spectral features of the disk emission, including the standard disk model, 
the inverse Compton scattering, and the slim disk model. In this paper, to understand the 
disk structure in the high state, we focus on whether the standard disk picture is valid or 
broken against the luminosity and disk temperature, by considering the effects of the inverse 
Compton scattering and the slim disk, rather than reclassifying the high states. Such an 
approach is important to describe many complex features found in the bright period of BHBs, 
not phenomelogically, but physically. 

3. Observation and Data Reduction 

The black-hole candidate 4U 1630-47 has been continuously monitored by the 
RXTE/ ASM (Bradt et al. 1993; Levine et al. 1996). Figure 1 shows time histories of a 
1.5-12 keV count rate and a hardness ratio of this source in 1996-2004, obtained with the 
RXTE/ ASM. The source showed five outbursts since 1996. We analyzed 322 pointing obser- 
vations with the RXTE/PCA of the five outbursts, which are indicated with arrows at the 
top in figure 1. In addition to the PCA data, the HEXTE data were analyzed only when 
their data quality was good. We analyzed the PCA standard-2 data and the HEXTE standard 
mode data for spectral fittings, and performed data selection on the condition of a minimum 
elevation angle of 10° above the Earth's limb. Since the spacecraft attitude is not stable just 
after changing the pointing direction, we excluded the data obtained when the actual pointing 
direction is more than O.T)2 away from the aimed direction. The selected data of the individual 
proportional counter units were co-added for spectral analyses. Systematic errors of 0.3% were 
added to the PCA spectra; furthermore, we added 2% and 10% systematic errors to the 4-8 
keV and 25-30 keV band, respectively, referring to the PCA calibration results in Jahoda et 
al. (2005). The latter is due to the uncertainty associated with the instrumental Xenon edge. 
The response matrices are generated by the ftool pcarsp, considering the change of calibration 
parameters, including the gain of PCA. 

We also utilized the ASCA/GIS data (Ohashi et al. 1996; Makishima et al. 1996) for 
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constraining the spectra at lower energies. The ASCA observation of this source was performed 
on 1998 February 26 (MJD = 50870), and this was almost simultaneous with the RXTE/PCA 
observation (observation ID 30188-02-20-00). The GIS events were extracted from a circular 
region of 6' radius centered on the image peak, after selecting good time intervals in a standard 
procedure. The dead-time fractions were determined to be 9% for GIS 2 and 12% for GIS 3 in 
reference to the count-rate monitor data (Makishima et al. 1996). 
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Fig. 1. The top panel shows time histories of the X-ray count rate and the hardness ratio (HR) (5-12 
keV/3-5 keV) of 4U 1630-47 by the RXTE/ ASM observation in 1996-2004. The PCA pointing observations 
are indicated with arrows at the top. "A", "G", and "ID 40112" indicate the observation data IDs 
referred to in the text. The bottom panels are enlargements of outbursts 1, 2 (left), and 5 (right). 



In pointing observations of ID 40112 in figure 1, the emission from 4U 1630-47 was very 
weak, and the count rate was constant, but higher, than the instrumental background level of 
the PCA. Considering further that a strong iron-line is present in the spectra, it is thought 
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that diffuse X-ray emission from the Galactic plane was mainly observed in this period. This 
emission must be a background for the observed PCA spectra of 4U 1630-47. We fitted 29 
pointing observations of ID 40112 with a thermal bremsstrahlung and a Gaussian representing 
the iron-line, multiplied by the interstellar absorption. Table 2 gives the values of the model 
parameters, each averaged over the 29 spectra. In a subsequent model fitting to the PCA spectra 
of 4U 1630-47, we included the Galactic X-ray background as a fixed model component, using 
the parameters in table 2. 



Table 2. Parameters for the Galactic diffuse background model, obtained from the data ID 40112. 



model wabs* 




bremss* 




Gaussian* 


iV H 


kT 


Normalization^ 


Line E 


Normalization 


(10 22 cm- 2 ) 


(keV) 




(keV) 


(photons cm -2 s _1 ) 


average 1.62 


17.0 


0.025 


6.58 


5.98 x 10" 4 



*: "wabs", "bremss", and "Gaussian" are spectral models in the xspec package; "wabs" represents the 
interstellar photoelectric absorption, "bremss" represents bremsstrahlung, and "Gaussian" represents an 
iron line. 



f: The normalization of "bremss" is in unit of 3.02 x 10 15 /(AttD 2 ) x J n e riidV [cm 5 ], where n c and m 
are electron and ion densities, V is the volume, and D is the distance to the object. 



4. Data Analysis and Results 

4-1- Spectral Fitting with the Standard Disk Model 

In order to grasp how the X-ray spectra of 4U 1630-47 varied with time, we fitted all 
of the spectra with the MCD model plus a powerlaw model, multiplied by the interstellar 
photoelectric absorption. We left all of the continuum parameters free to vary. Since the 
interstellar absorption toward 4U 1630-47 is not well known, we first fitted all of the spectra 
with the column density, Ah, left free, and then obtained an average value of Ah = 9.5 x 10 22 
cm -2 . We then repeated the fitting with Ah fixed to this value. The Ah values taken or obtained 
in previous studies are in the range of (5-12) xlO 22 cm -2 (Tomsick et al. 1998; Kuulkers et 
al. 1998; Cui et al. 2000; Dieters et al. 2000; Trudolyubov et al. 2001), and thus our 
value is consistent. The subsequent analysis employs the same value throughout. Because a 
broad iron-K absorption edge is often observed for black-hole binaries (Ebisawa et al. 1993), 
we multiplied the powerlaw component by a smeared-edge model (hereafter smedge model) to 
improve the fits around the iron-edge region. The edge energy, -E S medge, was typically found to 
be around 9 keV, and the optical depth, r sme dge, was 0.1-3, with rather large errors. The edge 
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width was not constrained. Therefore, hereafter, we do not discuss the edge feature. Hereafter, 
we call the above model as an "MCD plus powerlaw model" . 

The best-fit parameters for a typical spectrum observed in each outburst are given in 
table 3, where the above model is indicated as "MD" model. In figure 2, we show time histories 
of the X-ray luminosity, best-fit parameters, and reduced x 2 . Based on the identification of the 
soft/hard transition by Tomsick, Kaaret (2000), we here do not treat the data with an X-ray 
luminosity of < 6 x 10 37 erg s _1 , since such data are considered to be in the low/hard state. In 
addition, we discarded the data whose inner disk temperature is less than 0.8 keV, because the 
heavy interstellar absorption does not allow us to study such soft emission with confidence. In 
previous works, only outburst 2 was analyzed well, and our results concerning the MCD plus 
powerlaw model fits in the outburst 2 are almost consistent with those in Oosterbroek et al. 
(1998) and Tomsick, Kaaret (2000). 

The apparent inner disk radius, r- m , is calculated as Kmcd = 
[(r in /l km) x (10 kpc/D)] 2 cos#. The value of r in can be related to the true inner ra- 
dius, i?i n , as i?i n = ft 2 £r; n (Makishima et al. 2000). Here, k ~ 1.7-2.0 (Shimura, Takahara 
1995) is a spectral hardening factor, and £ = 0.41 (Kubota et al. 1998) is a correction factor 
for the inner boundary condition. We calculated the disk luminosity as L disk = 4:iirf n x aT^ 
from the best-fit apparent inner radius, ri n , and the inner disk temperature, T m . The powerlaw 
luminosity, L pow , in the 1-100 keV band is derived from the best-fit model, assuming an 
isotropic emission, and then the total luminosity is calculated as L tot = L disk + L pow . 

In figure 3, we show a scatter plot between Ldisk and T m in order to see the overall 
behavior of the X-ray spectra. The data points, defined as group X in figure 3, generally follow 
the relation Ldisk oc (the steeper solid line in the figure) which is expected for a standard 
accretion disk with a constant inner radius. However, two other groups of data deviate from 
this relation. One group (Y) of data shows Ldisk to be higher than 1.5 x 10 38 erg s -1 , and 
follow a flatter dependence on Ti n than the standard relation of Ldisk oc T^. The other group 
(Z) consists of those data points whose Ldisk is much lower than those of the former two groups 
for the same Ti n , or whose T\ n is unusually high for black hole binaries. Below, we examine the 
data of individual groups to understand their spectral states. 

Group X includes outburst 3 and the latter part of outbursts 2, 4, and 5 (hereafter 2b, 
4b, and 5c, respectively). These data points have a constant i?i n , though Ldisk changes from 
5 x 10 37 to 2 x 10 38 erg s -1 . The fits are acceptable, and as shown in figure 4c and table 3 
for outburst 2b, the spectrum is represented by a dominant soft component and a powerlaw 
tail. These results indicate that the standard disk picture is realized in these periods. We 
hence ascribe the group-X data to the standard regime, after Kubota, Makishima, Ebisawa 
(2001) and Kubota, Makishima (2004b). The true inner radius, i?i n , can be estimated to be 
30-45(D/10 kpc) cos# / ( l/v^) 2 km, corresponding to a non-spinning black-hole mass of 
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3.4-5.1(D/10 kpc) [cos0/(l/V3)J 2 M Q . These values are reasonable for a black hole binary. 

The former part of outburst 2 and the middle part of outburst 5 (hereafter, 2a and 
5b respectively) are included in group Z. Sometimes, the values of Ri n are less than 10 km 
and Ti n is also as high as 2 keV or more. These values are unusual for black hole binaries. 
The typical spectrum in outburst 2a is shown in figure 4b. The hard powerlaw component 
is dominant compared with those in the standard regime (figure 4c). The fits are relatively 
poor; the reduced x 2 is relatively large in outburst 2a and quite large, >2, in outburst 5b. The 
luminosity of the powerlaw component dominates that of the MCD component for almost all 
of the group Z data, as seen in figure 2. Outburst 2a contains the period of the very high state 
defined in Trudolyubov, Borozdin, Priedhorsky (2001). These properties make this outburst 
reminiscent of the anomalous regime of GRO J1655-40 and XTE J1550-564 (Kubota et al. 
2001, 2005). Hence, these data points are thought to correspond to the anomalous regime, 
which seems to appear at L tot > 2.5 x 10 38 erg s _1 for 4U 1630-47. 

Outburst 1 and the former part of outburst 4 (4a) are classified into group Y, so are 
two parts (outburst 5a) of outburst 5 before and after outburst 5b (figure 2). In outburst 1 
(figure 4a and table 3), the reduced x 2 value is significantly large, even though the spectral 
shape is represented by the dominant soft component with a weak powerlaw tail as well as the 
standard regime. Large reduced \ 2 values indicate that the spectra cannot be well represented 
by the standard disk model. As shown in figure 4a, the residual is significant between the data 
and MCD model around 5 keV. A narrow deep-edge structure is required around 7 keV for 
the powerlaw component in the MCD plus powerlaw model fitting, although this is thought 
to be artificial to explain the residual between the data and the MCD model, which is not 
the best model. Throughout this outburst, Ldisk is very high at > 1.5 x 10 38 erg s -1 , L pow is 
quite low compared with L tot (figure 2), and R in = 25 km is slightly smaller than that in the 
standard regime. Furthermore, the photon index of the powerlaw component is larger than 2.5, 
in contrast to the typical value of 2.0-2.5 in the standard regime of this source. These features 
of outburst 1 are the same as found in the apparently standard regime of XTE J 1550-564. 
Outbursts 4a and 5a show several features that are somewhat different from those of outburst 
1. The powerlaw component is not as weak as shown in figure 4d, and the photon index of 
the powerlaw component, 2.2-3.0, is only slightly larger than the typical values in the standard 
regime. Nevertheless, the disk luminosity Ldisk is very high at > 1.5 x 10 38 erg s _1 , and follows 
the relation of Ldisk oc T] 2 , as well as outburst 1. Therefore, while the spectral features of the 
powerlaw component are different among outbursts 1, 4a, and 5a, the behavior of the MCD 
component in these periods is similar to that of the apparently standard regime. 
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Table 3. The best-fit parameters. 



model 1 " N H * T in * R in * a ph ' N thc * L d i a k/itho/£ P ow * E smedge ' T amodgo ' W amodgo * X 2 /d.o.l 

(kcV) (km) r thc * (10 3s ergs _I ) (koV) (kcV) 

P* 



A (outburst 1) : 


apparently standard regime 
















MD 9.5 


, o 9 + 0.05 

1 - JZ -0.0l 


28.6+. / 7 


9 71 +0.29 
z ' ' 1 -0.74 


- 


2.24/-/0.11 


7 00+ 52 
' - uu -o.oo 


>1.00 


<0.01 


92.7/52 


PF 9.5 


1 36+°' 07 
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01-08-00(50899), 40418-01-03-00(51310), 60118-01-07-00(52001), 70417-01-02-01(52532), and 30188-02-20-00(50870). 
f:"MD" is a MCD model plus a powcrlaw model, multiplied by an interstellar absorption and a smeared edge (called as MCD plus 
powcrlaw model in the text). "CM" is a model in which the thcomp model is included in addition to the "MD" model. "PF" is a 
p-free disk model (see text subsection 4.3) plus a powerlaw model, multiplied by an interstellar absorption and a smeared edge. 
|:A^h is the column density of the Galactic interstellar absorption in unit of 10 22 cm -2 . Ti n and Ri n are an inner disk temperature 
and radius, respectively. a p h is a photon index of the powcrlaw component. TVthc and Tthc are a normalization and a spectral index 
of the thcomp model (see text subsection 4.2). p is a parameter of the p-frce disk model. I/disk, £ P ow, and Lthc are the luminosity 
of the MCD, powerlaw (1—100 keV), and thcomp (0.01-100 keV) models, respectively. £? S medge, T sme d ge and W S medgc arc an edge 
energy, depth and width, respectively, for the smeared edge model. 

§:Rcsults by fitting the PCA and HEXTE spectra. The data other than 'B' and 'G' are results by using the PCA spectrum. 
||:Results by simultaneously fitting the PCA, HEXTE, and ASCA-GIS spectra. 
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Fig. 2. Time histories of the spectral parameters obtained by fitting the PCA spectra with an MCD 
plus powerlaw model (MD model in table 3). Only the pointing observations of five outbursts are plot- 
ted. The light curves of outbursts 1 and 2 are shown in the top-left panel, those of outbursts 3 and 4 
are in the top-right panel, and those of outburst 5 are in the bottom-left panel. In each panel, time 
histories of luminosity L [10 38 erg s _1 ], Tj n [keV], i?; n [km], photon index a of the powerlaw, reduced 
X , and QPO frequency / [Hz] are plotted from top to bottom. In the time history of luminosity, green 
squares, blue triangles, and black circles indicate the disk luminosity Ldisk, the powerlaw luminosity 
L pow , and the total luminosity Ltot, respectively. The horizontal axis of the plot is defined as an offset 
day from the date that is specific for each outburst. The origin dates are MJD 50200, 50762, 51038, 
51420, and 51798 for outbursts 1, 2, 3, 4 and 5, respectively. For reducing the panel number, we plot- 
ted the outburst 1 and 2 in a single panel, and also did the outburst 3 and 4. The dashed ellipses and 
squares indicate the anomalous regime and apparently standard regime, respectively, while others are in 
the standard regime. For the labels "2a" , "2b" , "4a" , "4b" , "5a" , "5b" , and "5c" , see the text in detail. 
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Fig. 3. Ldisk plotted against T; n obtained by the MCD plus powcrlaw model (MD model). The blue, 
green, and red points represent group X. Y, and Z, respectively. Ldisk is in units of 10 38 erg s" 1 . Sec the 
text for the definition of the groups. The solid lines represent the Ldisk ^ ^in an d -^disk oc T? n relations, 
where the disk inner radius i?; n is assumed to be 46.1 km for the former relation. 
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Fig. 4. Examples of the spectra in different spectral 3tates. The data and their la errors are indicated by 
crosses. The solid line shows the best-fit model (MCD plus powerlaw model in the text, or MD model in 
table 3), which includes an MCD (dashed line), a powerlaw (dot-dashed line), and a Galactic background 
(dot line). Bottom graph of each panel shows residuals of the source counts from the best- fit model. 



4-2. Consideration of the Inverse Compton Scattering 

In the group-Z data, the spectral features are in good agreement with those in the 
anomalous regime of GRO J1655-40 and XTE J1550-564 (Kubota et al. 2001, 2005; Kobayashi 
et al. 2003). We therefore consider, after these previous studies, that the enhanced powerlaw 
component with a steep slope is produced when a certain fraction of photons from the optically- 
thick disk are Compton-up-scattered by hot electrons around the disk. Following Kubota, 
Makishima (2004b), we refitted these spectra with a three-component model, in which a 
thermal Comptonization component (thcomp; Zycki et al. 1999) is added to the original MCD 
plus powerlaw model. We hereafter refer this model to the "CM" model. We assumed that 
the seed photons are provided by the MCD emission, so the maximum color temperature of 
the seed photons was tied to T in . Furthermore, the electron temperature, T e , was fixed to 20 
keV. Hence, the thcomp model has two free parameters: a thcomp photon index, r tnc , which 
expresses the spectral shape below kT e , and a normalization, A^hc- Since the parameters of the 
thcomp and the powerlaw models can not be constrained independently, we fixed the photon 
index of the powerlaw component to 2.0, which is the average value in the standard regime. In 
figure 5, we show the time histories of the spectral parameters in the same manner as figure 2. 
Thus, the fit has been significantly improved, as shown in table 3 (see the line of "CM" model 
in table 3). 

Following Kubota, Makishima (2004), we plot L disk + L thc against T in in figure 6a, where 
L t h c is a luminosity of the thcomp model in the 0.01-100 keV. The values of T; n become 0.9-1.4 
keV, which are almost intermediate between the highest end of the standard regime and the 
lowest end of the apparently standard regime. However, the obtained values of Ldisk + L thc 
must be higher than the intrinsic luminosity of the underlying optically thick disk, Lf^ sk , due 
to the inverse Compton scattering. To correct this effect, following the equation by Kubota, 
Makishima (2004), we calculate the intrinsic disk flux as 

x ( p — ^) photons s -1 cmtl,) 

V 1 ke V / 

where F^ isk and F p hc are 0.01-100 keV photon flux of the direct disk component and the 
Comptonized component, respectively. Tj n refers to the disk temperature obtained by con- 
sidering the inverse Compton scattering. The intrinsic disk luminosity, Lf^, can be estimated 
as L^ sk = 4:iirf n x cxT^, where r in is given by equation (1). In figure 5, we plot a time history 
of -Ri n estimated from r- m obtained by this relation. Almost values of i?j n are obtained as 30- 
50 km, which are in good agreement with those in the standard regime. When the powerlaw 
photon index is fixed to a larger value of 2.3, the values of R[ n become smaller in such a way 
from ~40 km to ~32 km; the difference is not significant, and values obtained here are both 
close to that of the standard regime compared with the value in the MCD plus powerlaw model 
fitting. Figure 6b shows the estimated L^ sk against Tj n ; the anomalous regime data points now 
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align better along the standard disk locus, implying that the R m returns to its proper value. 
Thus, in the anomalous regime, the disk itself stays still in the standard state, whereas the 
inverse Compton-scattering process converts a significant fraction of its emission into the hard 
component, as suggested in GRO J1655-40 and XTE J1550-564. 

Even after these corrections, some data points in figure 6b deviate from the standard 
relation. They are from outburst 5b, where the total luminosity becomes the highest among 
the RXTE observations. We speculate that the underlying disk, itself, becomes deviated in 
these period from the standard picture. 

In outburst 2a, we have a simultaneous ASCA GIS coverage at MJD = 50870. The 
GIS covers the lower energy range of 0.7-10 keV than the PCA, which covers the energy range 
of 2-60 keV. Therefore, we can obtain more accurate values of T m and A% by using the GIS 
data. We included the 1-10 keV GIS spectrum in addition to the PCA and HEXTE spectra, 
and fit them simultaneously with the above "CM" model, without fixing the value of Ah- The 
results are given in table 3 (label "G"), and the achieved fit is shown in figure 7. In spite of 
the wide energy range (1-100 keV) and the extremely high data statistics (of the GIS and the 
PCA/HEXTE), the "CM" model has provided an amazingly successful (x 2 /d.o.f = 191.2/210) 
joint fit to the spectra. The best-fit parameters have remained close to those obtained with 
the RXTE data alone, without any knowledge for energies below ~ 3 keV. These results give 
a strong justification to the "CM" modeling. As can be seen from figure 7a, the flux in an 
intermediate energy range (5-10 keV) is now carried by the thcomp component, with the MCD 
and powerlaw components dominating in the lower and higher energies, respectively. This makes 
a contrast to the conventional MCD plus powerlaw fit to the same PCA spectrum (figure 7b), 
in which the model tried to reproduce the intermediate-range flux by increasing T m (which in 
turn reduced i? in ) and steepening the powerlaw. 
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Fig. 5. Same as figure 2 concerning the standard regime and the apparently standard regime data, but 
results in the anomalous regime are replaced by those obtained with the CM model. In each panel, time 
histories of luminosity L [10 38 erg s _1 ], T; n [kcV], R[ n [km], spectral index r t hc of the thcomp, fraction of 
ithc against Ldisk + ithc, and reduced x 2 arc plotted from top to bottom. The horizontal axis represents 
the same as that in figure 2. 
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Fig. 6. Same as figure 3, but the data points in the anomalous regime are re-calculated as (a) idisk + ithc 
or (b) i d n i* k , where L 1 ^^ is estimated by considering the conservation of the photon number. The longi- 
tudinal axis is in unit of 10 38 erg s _1 . 
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Fig. 7. (a) Joint fit to the GIS, PCA, and HEXTE spectra obtained simultaneously on MJD=50870, 
with the "CM" model. For the PCA spectra, the best-fit model and individual three additive components 
(MCD, thcomp, and powerlaw) are plotted. Only the best-fit model is shown for the GIS and HEXTE 
spectra, (b) Joint fit to the PCA and HEXTE spectrum on MJD=50870, with the MCD plus powerlaw 
model. Here, for simplicity, only the PCA spectrum is shown. The solid, dashed, dot-dashed, and dot 
lines, together with cross points, are the same as figure 4. 
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Here, for simplicity, only the PCA spectrum is shown. 
4-3. Evidence of Slim Disk 

The spectral features in outbursts 1, 4a, and 5a (group Y) apparently resemble those in 
the standard regime. However, these data points follow the relation L disk oc T 2 n , as noticed in 
subsection 3.1, and the spectral shape is not necessarily represented by the standard MCD plus 
powerlaw model. This suggests that the spectral state in group Y corresponds to the apparently 
standard regime (Kubota, Makishima 2004b;Kubota et al. 2005), in which the inner disk may 
be described by the slim-disk solution. In order to examine whether outburst 1 data can be 
described by the slim-disk picture, we here employ the p-free disk model, which mathematically 
extends the MCD model (Mineshige et al. 1994; Hirano et al. 1995; Kubota et al. 2005). 
This model assumes the radial temperature profile to be as T(r) = T in (r/r in )~ p , where p is a 
free parameter. The profile with p = 0.75 corresponds to the MCD model, while p <0.75 yields 
a spectrum softer than the MCD model for the same Tj n . That is, a small p suggests a low 
radiative efficiency at the inner portion of the accretion disk. 

Following Kubota, Makishima (2004b), we fitted the spectra of the apparently standard 
regime, by replacing the MCD component with the p-free disk model, by fixing the powerlaw 
photon index to 2. We hereafter refer this model to the "PF" model. The free parameters 
are T- m , p, normalization of the p-free disk and the powerlaw model, and the smeared-edge 
model parameters. Here, we analyze the data of outbursts 1 and 4b; the former is in the 
apparently standard regime, while the latter is in the standard regime for comparison. The 
best-fit values and the reduced \ 2 are plotted in figure 8a. It can be seen that the reduced 
X 2 for outburst 1 (apparently standard regime) dramatically becomes much smaller from ~ 4 
to ~ 1.5 (see also figure 2 and table 3). The reduced \ 2 values are still somewhat large, 
probably due to the incompleteness of the p-free model, which does not take into account the 
inner boundary condition and relativistic effects. For outburst 4b (standard regime), the R m 
values widely scatter with large errors, indicating that the p-free disk model introduces an extra 
parameter. Figure 8b shows the best-fit values of p, plotted against T in obtained with the MCD 
plus powerlaw model fit. In this diagram, the values of p in the apparently standard regime 
are 0.5-0.7, which is smaller than ~0.8 in the standard regime at T- m ~ 1.1 keV. Therefore, 
we suggest that the radiative efficiency at the inner accretion disk is lower in the apparently 
standard regime than in the standard regime; another cooling process, like advective cooling, is 
required in addition to radiative cooling. The theoretical solution of the slim disk takes into 
account the effect of advective cooling, and predicts smaller values of p than the standard disk 
(e.g., Watarai et al. 2000). In addition, a relation of r in oc TJ" 1 is found, when the spectra 
predicted by the slim disk solution are fitted with the MCD plus powerlaw model (Watarai 
et al. 2000). This explains the relation of Ldisk oc T? n of group Y. Thus, the X-ray spectra 
in group Y show evidence of a slim disk. This is the second example of such a state after 
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XTE J 1550-564. We note that the value of p at a lower T; n in outburst 4b (standard regime) 
is somewhat smaller, and such a trend is reported to be due to because the MCD model is a 
mere approximation of the exact standard-disk solution (Kubota, Makishima 2004b; Kubota 
et al. 2005). 
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Fig. 8. (a)Time histories of the best-fit parameters and reduced x 2 for the "PF" model fits. The top and 
bottom panels are for outbursts 1 and 4b, respectively. (b)Bcst-fit values of p in the p-free disk model, 
plotted against Ti n obtained through a spectral fitting with the MCD plus powerlaw model. The circles 
and diamonds correspond to the apparently standard regime and standard regime, respectively. 
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5. Discussion 



5.1. Three Spectral States Concerned with the Disk Structure 

Through an extensive spectral analysis of the RXTE data covering the five outbursts in 
1996-2004, we confirmed that most of the X-ray spectra in the high state of 4U 1630-47 can 
be classified into three distinct spectral states. One state is well explained by the standard 
accretion disk, which appears when T m is less than 1.2 keV, or the X-ray luminosity is below 
2.5 x 10 38 erg s _1 . In this standard regime, the data satisfies the Ldisk oc relation, and thus 
this regime corresponds to the standard high/soft state. In this section, we give answers to 
aims 1, 3, and 4 of section 1. 

When the X-ray luminosity exceeds 2.5 x 10 38 erg s -1 , the other two states appear. One 
of them, the anomalous regime, is characterized by the dominant steep (T ~ 2.5) powerlaw, as 
well as the unusual values of R m and Tj n , which are obtained when fitting the spectra with the 
MCD plus powerlaw model. The dominant powerlaw component can be understood in such 
a way that a significant fraction of the disk emission (soft component) is converted into the 
hard component through inverse Compton scattering. The apparent variability of i?i n can be 
explained by this effect, and i?; n remains constant when we consider the Comptonized photons 
for the disk flux. Thus, after GRO J1655-40 and XTE J1550-564, 4U 1630-47 is the third source 
to exhibit the anomalous regime. The data of this regime is mostly in the very high state or 
intermediate state. Therefore, our interpretation by the Comptonization picture indicates that, 
in some parts of the very high state or intermediate state, the disk structure is still standard, 
but the ambient hot electron corona that Comptonizes the disk emission is formed and evolved. 
The quantification of spectral structures in the framework of Kubota, Makishima, Ebisawa 
(2001) is important, since the black hole mass measurement is still available, even in the very 
high state, and we can constrain the solid angle of the hot electron corona around the disk from 
the fraction of Comptonized photons in the disk emission. 

The other state appears in outbursts 1, 4a, and 5a, where the spectral features are 
apparently similar to those in the standard regime. However, the disk luminosity is found 
to follow the relation L^k oc T?, indicating a low radiation efficiency of the accretion disk. 
This feature is similar to those of the apparently standard regime observed in XTE J1550-564 
(Kubota, Makishima 2004b), and maybe interpreted as the formation of an optically thick 
advective disk, namely, a slim-disk. The apparently standard regime also corresponds to the 
standard high/soft state with very high X-ray luminosity. However, X-ray spectra are not 
exactly represented by the standard MCD plus powerlaw model, and prefer the p-free disk 
model; and thus, we here stress that the disk structure somewhat deviates from the standard 
disk. 

As above, we can mostly understand the disk structure in the high state of 4U 1630- 
47 within the framework of Kubota, Makishima, Ebisawa (2001) and Kubota, Makishima 
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(2004b), after GRO J1655-40 and XTE J1550-564. However, some parts of the spectra cannot 
be well explained by the above three spectral states, and thus further studies are needed to 
understand the disk structure of a very high state or an intermediate state completely. 

5.2. Flux Fraction of the Powerlaw Component and Critical Luminosity 

Although we have utilized luminosity-dependent spectral changes to distinguish the three 
states, there can be a simpler way of classification, which is usable even in limited snap-shot 
observations. In figure 9, we plot the luminosity ratio of the powerlaw to the total component, 
L pow /L tot , against the total luminosity, L tot , where L pow is that obtained in the MCD plus 
powerlaw model fitting. It can be seen that the three states are clearly separated in this 
diagram. The standard regime (blue points) locates at L tot < 2.5 x 10 38 erg s _1 . The anomalous 
regime (green points) is characterized by L tot > 2.5 x 10 38 erg s _1 and L pow /L tot > 0.5. The region 
of L tot > 2.5 x 10 38 erg s _1 and L pow /L tot < 0.5 corresponds to the apparently standard regime 
(red points). The threshold luminosity of 2.5 x 10 38 erg s _1 should be different among black 
hole binaries, and hence cannot be constrained from the limited data. However, the anomalous 
regime can be identified only by L pow /L tot . Therefore, if the data include the anomalous regime, 
we can constrain the threshold luminosity, and thus identify the apparently standard regime. 
This suggestion can be tested by using the RXTE data of a number of black hole binaries. 

Throughout the five outbursts of 4U 1630-47, we have found that the well-defined critical 
luminosity of L crit ~ 2.5 x 10 38 erg s _1 divides the standard regime from the other two, with 
a good reproducibility. This is the first case where the critical luminosity L cr i t is almost the 
same among outbursts of one object, indicating that L crit is proper to the object. Therefore, 
objective 2 in section 1 is achieved as above. We may also estimate the Eddington luminosity 
as Le ~ (6.8-10.2) x 10 38 erg s -1 , using the black hole mass estimate of 3.4-5.1 M & (subsection 
4.1). Then, we find L cvlt /L E = 0.25-0.35 for 4U 1630-47. This value is somewhat higher than 
the values of L crit /L E ~ 0.15-0.20 observed from GRO J1655-40 (Kubota et al. 2001) and 
XTE J1550-564 (Kubota, Makishima 2004b), while lower than that for LMC X-3, L crit /L E > 
0.9 (Kubota et al. 2001). The critical temperature corresponding to L crit is ~1.2 keV for 4U 
1630-47 and GRO J1655-40, and 1.0 keV for XTE J1550-564. Therefore, the transition is not 
determined only by the temperature. We speculate that L crit and the transition temperature 
may be determined by the combination of the black hole mass, spin, and the accretion rate. 

5.3. Issues on the anomalous regime 

Since the anomalous regime is important to understand the unified view of the very 
high state or intermediate state, we here briefly discuss two unresolved issues concerning the 
anomalous regime. 

Outburst 2 was accompanied by an optically thin radio flare (Hjellming et al. 1999). The 
radio flare lasted from MJD 50840 to MJD 50880, in good agreement with the anomalous regime 
period (figure 2). XTE J1550-564 also exhibits an optically thin radio flare (Wu et al. 2002) 

21 




L tot (10 38 erg s X ) 

Fig. 9. Luminosity ratio of the powcrlaw to the total component obtained with the MCD plus powcrlaw 
fit, plotted against the total luminosity. The luminosities are in units of 10 38 erg s _1 , and refer to values 
in figure 2. 



in the anomalous regime (Kubota, Done 2004a). In addition, for GRO J1655-40, XTE J1550- 
564, and 4U 1630-47, the radio flare and QPOs are associated with only the anomalous regime. 
These results suggest that the optically thin radio flare and QPOs have relations with the 
Comptonizing plasma, rather than the accretion disk, and these two phenomena are thus good 
indicators of the anomalous regime. Further studies of the correlation between the anomalous 
regime and a radio flare are interesting from the view point of the disk/jet connection. 

In figure 6b, most data points of outburst 5b, which we identified as the anomalous 
regime, deviate from other data points toward higher luminosity. Such a deviation is similar to 
the case of the strong very high state in XTE J1550-564 (Kubota, Makishima 2004b). They 
indicated that the accretion disk is not in the standard state and/or does not reach the last- 
stable orbit. In the case of 4U 1630-47, the high state continued just before outburst 5b, and 
thus the situation is different; the truncation of the disk is possibly caused by extremely high 
radiation pressure. Since the total luminosity was the highest in this outburst, the radiative 
pressure became highest. Then, the inner disk material is thought to produce a hot and optically 
thin corona. As a result, the disk is truncated and the inner disk radius, R{ n , becomes large, 
leading to the higher luminosity. 
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